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ABSTilf\CT 

The   stai  ility   characteristics   of   several   .  aroLropic  inouels 
(the   oivergent,    one-parameter  and    the   stratified  models)   are    inves- 
tiijated    foi.    vari..'UQ    .>>uimctrical   .■■jmxi    ./tad   fields,    iucluui.iij,    :iaj,:.c- 
jet  and    aoulie- j  ;..t   profiles.      Two   differeat   approximations   are   used 
for    the   amplitude   of    tlie   per L'ur- atioiiL   stream  fuuction,    one    of    .  hich 
i'     .;asec   oa  a  i-i^lhou   oi    li-a.-nj   differences,    and   tr;u   secotiv;   on  a 
represeutatioii  '/y  a   finite   Fourier   serius.      The  results   of   the   t.vo 
methods   are  compared   to  one  anotiier  as    ;ell   as    to   those   of   previous 
.i..:i\ estimators  . 

The   second -o;..'jer   effects   of   the   unstable   perturi-ations   Oii  a 
doul'ie-jet   zonal    flo."  are  determined   from  tiie   converi^ence   of  meri- 
dional   transport    of   laoiuentuui   i^y    the   eddies.       It    '..:    sho:/a    tiiat    t.;it^ 
type   of    ;:onai    .viad    profile   chaji^es    to   a   Gin^le-jet   mean   f  lo  ;    iii  a 
fe\\?  i^ays  . 
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1.   latroductioii . 

The  iarotropic  model,  conceived   y  Rossi^y  raore  than  I  wo 
decavles  a^o,  tias  pla^eu  a  roie  in  aui'iitsricai  I'orecastiag  far  beyoiicl 
expcctatioiuj ,   iJheu  operational  forecasting  -.-'ith  the  Larotroplc 
model  began  in  tVie  early  'SO's,  it  .as  expected  chat  the  more  so- 
phisticateu  aroclinic  types  would  soon  replace  the  former;  however, 
difficulties  \vith  these  models,  to^-icther  .vith  furth.er  improve- 
ments in  the  barotropic  model,  have  caused  cootiaued  interest  ia 
the  latter. 

One  aspect  of  i^arctropic  flo..^  v/hich  has  been  investi^^ated  in 
some  detail  is  the  dynamic  stability  of  ::;onal  currents,  i,u.-a>  ly 
by  Kuo  [i]  and  Eliasei>  [2]  for  non-divergent  flo./,  and  for  the 
diver^^ent  case  by  Wlia—Mielsen  [j]  , 

la  the  first  part  of  this  paper,  the  stability  propertict>  ul 
varioc^s  larotroplc  models  v.'ill  be  studied  by  t;;o  different  ap- 
proaches, namely  a  finice-dif f erence  metiiod  and  a  finite  Fourier- 
series  method.   A  varxuty  of  :-.ouai  ■.;ind  profiles  'irill  Le  con- 
sidered, includin;-j  single  and  douole  wind-speed  maxima.   The 
results  of  the  t/o  methous  are  corapa'red  to  one  another  as  sjel'i   as 
to  those  of  previous  inves Liu^ators  . 

The  performance  of  the  numerical  predictio.i  models,  ■..hen 
integrated  over  periods  of  time  up  to  several  days,  .-eeks ,  or 
e  en  more,  has  been  recently  of  special  interest,   Thompson  \_-i} 
has  developeu  a  heuristic  tb^eory  of  long-period  velocity  varia- 
tions in  barotropic  flo;.   I,i  particular,  it  wab  sho-.m  from  this 
theory  that  a  double-jet  maxima  in  the  iional  Vv'inu  profile 


develops  in  t^o  or  three  days  from  an  initial  single- jet.   Such 
conditions  seem  to  be  ciiEracturis t ic  w  :  certaia  t}pcs  of  atmospheric 
currents  ./heu  a   loci;  occurs.   This  theory  is  Ijasec  upon  the  in- 
teraction bet  vcen  eociies  and  the  mean  flo/. 

Wiin-Nielsen  [3]  has  recently  sho  .'li  that  the  splitting  of 
a  sin,j,ie  isotach  maximum  in  the  meridional  v;inu  profile  (sia,.ile 
jet)  into  tv;o  ma.iinia  (double  jet)  is  a  dynamical  couseq-ence  of 
the  convergence  oi'  iiieriuioaa  I  laouientum  traaspoi't  .-y  tnt;  euLiiiJb, 
provided  there  is  barotropic  instability  in  an  initial  flo.v.   It 
may  he   noted  also  that  the  time  scale  of  the  chaii;^es  in  the  mean 
;;onal  flo..'  is  the  same  in  uOtli  tlieories,  as  may  ;.e  expecteu,  since 
the  physical  model  is  basically  the  same. 

Tlie  second  part  of  tuis  paper  is  devoted  to  a  study  of 
second-order  effects  of  the  perturi:ations  on  mean  ?.onal  flow 
which  has  double  maicima . 


2.        Stability   properties   of    the  diver  j,ifc~parameter  mooe  . . 

A.    jiasic   eq'-.ations. 

Tlie  vorci.clL_    eq^i.  or    the   uivergunt,    one -parameter  mouui 

(or   equivalent   baj.ctropic  model)  may  be    /ritten   in   the   foliowin;^ 

form  di.e  to  vJii.i--Mielse;i  j^..  ]  : 

|-{vY}+v-v(vV+f)=  q'H^ 

where  If^'  is  the  stream  function  j 

%   =  IK.   '-  VV^    ,  the  horizontal  uo. 

and  q^=  -i  4^  iJi_-  .  :) 

In  equation  (2.2)   (7=— f^)^^^    i-'  a  measLire  of  static 

stability,  Jq     a  standard  value  of  the  Corioil.s  j^ara.;:. ,  . 

A(p)  is  an  erapirical  function  describing  the  vertical  varia 

the  horizontal  v;ind  near  the  level  \;here  the      ,      ^  :.  _   -:  is 

applied  (500  mb),  namely 

V  (X,  y,  p,  t)  =A(p)V3  (X,  y,  t).  (2.3) 

By  applying  the  perturbation  m(Sthod,  the  prognostic  eqv^ation 
(2.1)  may  '  e  'Written  in  the  form 

BTir  , 

v;here  [J(\J)=  —    — ^  is  the  mean  i:oual  wind  velocity,  anu  "ij/" 

is  the  perturbation  stream  function. 

The  perturbation  stream  function  ■..i.j.x  :.e  contiicered  in  the  form 


y\x,  '-j.t)  =  <i)(y}^'>^^^-^^> 


(:..  ) 


SuL^stiLuting  {2. 5)    into    (2.4)  yields  the  result 


(0-c)[^-./^l,[e-g^fc-Jt=^o   . 


(2.6) 


27T 
la  (2.5)  M  =  ~r'"  (^  ^^'^^  quantity)  auu  c,  the  phase 

velocity,  may  .'e  complex,  that  is 

C   =  Cr    -\-    i  C^        ^       .  (2.7) 

If,  in  (2.7),  Ci  ^  O     the  amplitude  v/ill  ■■e  an  exponential 
fuuction  of  time   t  aaC   the  ',/a\'e  is  saiu  to  'e  unstal'le. 
For  dye       >    ^^'^'^^   arapiitude  will  increase  -Lth  time;   if,  on 
the  other  iiau-. ,  Ci   -C  Q      ,  the  .,ave  will  be  damped, 

Equatlo.-;  (2,u)  ia  a  second-order  differential  equation  an,,  ^l. 
to   e  solved  for  certain  ivoundary  conditio      .-ince  disturbances 
in  :ionai  curreiits  are  normally  bounded  to  the  north  and  south,  the 
simplest  type  of  suitaMe  ioundary  cono^u. oi.,   ._,   ..:   rotten  Li. 
the  form 

(f)(0)=  (j)(D)  =  0  ,  (2.8) 

corresponoin;^  to  ri^i^id  boundaries  at  y  =  o  and  y  =  D,  where  the 
perturbation  stream  fuiiction  is  assumed  to  vanish. 

No-;,  as  au  >:.i.;5  en -value  problem,  we  are  interested  i...  ^.u- 
relationship  betv/een  c,  jj.     and  other  pertinent  parameters  for 
v/hich  there  exists  solutions  of  the  system  (2.6)  and  (2.8). 

Since  t.ic  ut'.^.lc:  current  is  a  function  of  y,  the  mathematical 
problem  here  is  rather  complicated;  and  generally  v/e  cannot  find, 
even  for  simple  functions  U  (y),  an  analytical  expression  for  the 
pertur:atLon  amplituce  p  (y).   b.ence  it  appears  logical  to  atte.upt 


solution   1  y  a   fiuilic-olf ference  metho-.;   or  a    fiaite   Fourier- 
series   method    ./hicii  ca  i      t-  aciapte..    Lj    uhe   particular  loundary 
coaditiotis    (2.3), 


B.   Method  of  finite  (■inferences,  • 

In  the  atrnosphe...^,  .'^^locity  profxteo  u  ■_  li.u  .esterlies 
Seaerally  have  apprcciaiJie  symmetry  ./itb  respect  to  some  latitude. 
For  simplicity,  particular  symmatric  profiles,  especially  single- 
jet  and  oou-  j.e-_,'^L.  ^yrofiles,  './ill  be  discussec.. 

Since  only  symmetric  wind  profiles  and  >  ounuary  conditions 
are  considered,  sytTinietry  of  the  amplitude  0(y)  of  the  perturba- 
tion stream  f u;.L.i-i.o-i  i.iay  i:u  expected.   The^-tj..:jre,  only  half  oi 
the  v.'ind  belt  aee;.  i,e  considered. 

la  so  tioing  .'e  shall  consider  the  fuactiOiL  p(y)  as  ^iven 
i.y  its  values  ^^.    aL   i\   pc}jju.s  in  the  internal,   Q^V-O       > 
such  that  the  iiitoival  is  divided  by  these  points  into  N  +  1 

equal  parts,  e^ch  of  tlie  lenr,th  d  =   D    (figure  ]  "• .   The 

N+1 
meriaioaai  uerivatives  of  p(y)  v;iii  ue  approximateu  ;;y  the  cor- 
responding quotients  of  centered  differences  such  as 


fSUP) 


dy^^n  ^2  ~     .  .  (2.:) 

In  ^>eaerai  it  vnay  ..e  expected  tiiat  tiie  approximation  of  tiie 

derivatives  l.y  tliis  difference  ratio  ./ill  give  better  results 

as  the  value  of  N   increases,  althou;sh  the  mathematical  analysis 

becomes  more  tedious.   Ho.'ever,  for  simplicity  iii  this  sectioii 

only,  the  cases  of  N  =  5   and  N  =  7   will  l.e  considered. 

Fart  I.    K  =  5,  (or   .    -^  ) . 

o 

By  combining  (2.9)  and  the  given  boundary  conditions  (2.8), 
the  derive  lives  of  the  pertviriatiou  amplitude  at  each  point 
(':'.g.  1)  can  '^e  ■.written  in  the  form 


2t, 

-2to 

4>e  + 

•4^.-2 

t. 

t/    - 

d" 

(2.10) 


Since  (2,6)  iTiust  be  satisfiea  at  eacn  point  of  fig.  1,  the 
follov;in3  set  o£  equations  with  the  use  of  (2.10)  are  obtained: 

[(Uo-c)  (-z-im  ^d"(^-v:')  +  :j"q'c]<t>,  +  2  (Uo-c)ct>,     ^  0 

III)  1 :  : ! 

'.-.'here     U    ,    U.  j    U„ ,    U      ,    U,      and  IL,      are   the    v/alues   of     U     and  U 
o         JL         «:-         G  i  ^ 

at  V    5      y,      and  >%, ,      respectively.      Here   tlie   douliie  prime   repre- 

■"0     1       Z 

sents  the  second  derivative  with  respect  to   y.   For  a  known 

initial  rional  v/ind  profile,  (2,11)  constitutes  a  system  o'  ! 'Tee 

equations  in  three  unknowns  0  ,  0^  ,  and  p„ .   Since  the  systeai 

o   1       2 

(2.11)  is  homogeneous,  the  possiMe  values  of   c  rnay  iiov/  be 
obtained  from  the  co-.idition  that  the  determinant  of  the  syste.a 
(2.11)  must  vanish  in  order  to  have  non-trivial  solutions.   This 
conditions  leads  to  the  so-called  frequency  equation  in  the  form: 


cU(A-3)]^cU{ir>^iiU^(VMj,)(i-y^)] 
-+(y^-'-/)(a-Uor)+2v?fu>Li,)-Z'(/i-ur)]+c[-^{Tj,Dz(K' 

-/)4-T^"R-r(ll[L-hBU,)K{^(i^-^i^)+fD.+il2)r/-r/i)l        (^-12) 

{Vo-u.r}  -2kdoO,fZCu.TU,nrz-iJi.Y')]  +  [{T},uz{r-n 

■/here  ,  -,  „ 

1^,=  d^(:?-TTo)         , 


dV^-n,") 


and  '^=ry+jq-  ^ 

It  should  be  iiotec.  that  0     and  <5   ..ere  a^ouiue-  constant  iri 
deriving  (2.12).   Since  tlie  frequency  equatio.i  (2.12)  is  a  culic 
equation  in  the  phase  velocity  c,  in  general  there  existf^  three 
roots,  ;.'hich  may  iiicra<ie  a  pair  of  complex  roots.   The  ^i.:^  ^  i.iiai_> 
part  of  the  phase  velocity,  c.,  v/nich  is  necessary  for  unstable 
perturbations,  v;ill  be  considered  of  primary  interest  in  this 
paper. 

The  frequency  equation  (2.12)  has  been  i;olved  numerically 
for  a  number  of  cases  ..'ith  a  Cotitrol-Data-Corporation  1604 
computer,  using  tne  so-called  Huller  metlioa  \^C\   for  finding  the 
roots  of  polynomials.   For  each  of  the  initial  ./ind  profiles 

and  meridional  distar^ces  betweeii  the  rigid  '..boundaries,  the  solu- 
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tion  v.'as  obtained  for  various  ional  waveiengtns  (  1  rr—   ), 


keeping  ' ^'     aad  q   coustaut.   Then  the  effects  of  i^       and 
divergence   q   ./ere  teste  '  i  y  making  each  of  these  parameters  /.ero 
ill  a  vSeries  of  conii^utatioaK . 

It  is  of  special  interest  to  consider  the  stability  problem 
for  those  types  of  initial  wind  profiles  ./hich  are  not  only  sym- 
laetric  v/ith  respect  to  the  center  of  the  channel,  ',ut  also  fre- 
quently found  in  the  actual  atmosphere. 
Example  1.   Single-jet  ^asic  current. 

For  the  first  example,  the  :-onai  wind  profile  v/ill  be  con- 
siuered  in  the  forra 

U=  B  C/-fOS^'^')  (2.13) 

'..'hich  is  shovm  as  curve  (1)  in  figure  2.   Inserting  (2.13)  into 
(2,12)  and  varyiug  the  pertinent  parameters   B,   D,    A'    ,  a.;. 
q    leads  to  the  results  j^iven  in  figures  3,  ■'. ,  ^,    u,  ana  7. 

In  figure  3  the  meridional  width  of  the  chianuel  ^  D,  is 
2000  km  and  the  rnaximuvn  oonal  speed  at  the  center  of  the  channel 
is  60  m  sec    (E  =  30  m  sec   ).   Thi-,  case  corrcsponus  to  a 
rather  intense  single- jet  ./ind  profile,  ..'liich  is  often  found  at 
the  500-mb  level.   T'ne  ai-scissa  in  the  figure  is  conal  wavelength 
in  units  of  tiiousands  of  km,  v^yhilo  the  ordinate  is  the  imaginary 
part  of  the  phase  velocity   c.  .   The  figure  cc.ttains  five  curves 
corresponding  to  the  cases:   curve  (i)  ^  ~  0,    q  -  0,  (2)  ^/^  =  ^i^^   , 

q  =  0;   (3)  ,^=-/^4.,r   ,  q^=  0.i>  x  lo"^^  m  "'^ ;  (4)  /9  =  ^.^ 

''  -12-2  2  -12   -2 

q""=  1.0  X  10   M    ;   and  (5)  3  =  p^^   ,  q  =  1 . 3  x  10   m    . 

Here   /^^    represents  the  meridional  derivative  of  tliC  Coriolis 

parameter  at  latit.^de  45N   ^^^hich  has  the  value  1.618x10   m  sec 


By  compiirins  the  two  first  curvea,  it  i."  seen  that  the  para- 
meter ,6'  luis  a  3tabilir:iii5  Influeucii  on  the  pertur'ations  ia   the 
seiise  that  tlie  ina'-niitu'-ie  of  c,   is  svnalier  auu  the  hand   of  wavc- 
lep.gths  for  whicr.  iastability  occurs  is  more  uarro  /  for  B  ^  0 
tb.an  in  the  case  of  S    =   o.   Moreover,  !,y  compa--..,  -.le  first  t.  o 
c^-rves  jLth  th.e  lact  three  carves  it  is  also  seen  that  q,   which 
ir-:  a  r^ieasure  of  the  intensity  of  the  divergence  implied  by  this 
raodel,  h.as  further  stai..ili:;ing  effects.   Therefore,  the  divergent, 
one-parameter  model  is  more  staf'ie  than  the  non-divergent  model. 
Furtherraore ,  th^e  stabili:-:;ing  effects  of  Q   a.id  the  divergence  term 
are  much,  more  pro.iouacod  for  relatively  !>;.  ..■  .avelengths  in  the 
previously  unstable  bancs.   Thus  the  riiaxinium  c.'s  occur  at  shorter 
"..'ave lengths  f or  A' :^  o  '•i"'':''  'i i=  0  ,  than  in  the  case  of  ^  =  q  =  o. 

It  should  l^e  noted  lit  this  figure  that  the  instability  occurs 
only  for  the  intermediate  ;:;onal  \vavelength3  (in  this  case  2000-GOOO  hm 
for  A   =  q  =  o),  vhich  is  in  agreement  v/ith  the  average  lengtii  of 
observed  synoptic  waves.   In  other  .'ords,  for  sufficiently  short 
and  long  wavelengths,  instability  does  not  exist; .   This  is  similar 
to  the  res  .Its  of  the  t\;o-level  baroelinic  instability  discussed 
by  Thompson   [z]  . 

The  preceding  results  are  dicely  verified  (5y  figures  4,  5, 
and  6,  whicl\  contain  similar  computations  .7ith  the  meridional 
widths,  3000,  4000,  anu  6000  km,  respectively.   The  additional 
information  shcv;n  in  tliese  figures  is  that  the  band  of  unstable 
'..'ave lengths  is  shifted  to'.'ard  longer  wavelengths  such  that  c. 
becomes  maximum  at    i  2.0  for   6     =  <\    -  o\    and  also  the 
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magnitude  of  the  maxinrar.i  c.   tloes  aot  change  '.'ith  D  in  this  cast. 
Thio  result  is  coraplelely  cunsi  f^touL   Itii  cliasc.  [ ']  .   iloL'fcover, 
it  may  i^e  note.,;  that  the  atal^iii^iii^^  effects  of  3   aac  q  are 
generally  more  pro-aounceci  with  increasing  D;   tlras  the  >/aves  are 
more  stable  in  the  cliverr^ent  model  with  incrGf.r- '  .1-  '/idth  or  the 

cl.anaei.   For  example,  ';/ith  D  =  4000  i:ra,  all- -vaveleniiths  are  sta  le 

2  -12   -2 

in  the  case  oi:  R  —  ^^     and  q   =  1.5  x  10    \:\      ;    and  only  the 

..■aves  in  a  narrow  band  around  5000  kra  are  uristal-l'^  :^0'  ^  =  Btt.^ 

2  -17   -2 

and  q   =  1.0  x  10    vr;    .   For   D  =  6000  i;tn,  o;.iy  the  case, 

S   =  0     and  q  =  0,  is  unstable. 

Figure  7  contains  five  curves  giving  th-  "'  -^f.'iinary  root   c. 

as  a  function  of  t'lie  meridional  v;idth  ,D,  for  a   fix(=d  value  of 

the  2onal  v;avelengti;,   L  -   5000  km. 

Example  2. 

For   tlie   secon."   profile,    it  v/ill   be  assumed    that 

U     =     2B 
o 

U^   =       B  (2.14) 

u     =       " 

I 

•..'here     B   =  30  m  sec 

Usin^  the  finite-difference  approximation  for  the  derivatives  of 

?:onal  v/ind  speed  uith  respect  to  y  yields 

n:'=  0 
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lusertiu^  (2.14)  anu  (2,15)  into  (2.12)  a. id  performiag  the  accessary 
computations  leads  to  the  result;.  ,:..,.   ^  fij^ures  8,  9,  air;  10, 
corrcspondlnii  to  D  -  2000,  4000,  and  5000  I'm,  i-espectively . 
It  ij  readily  seen  that  the  results  of  the  instability 
faatLres  for  this  profile  are  quite  similar  to  those  of  the  first 
proJrile,  .'hich  ij  to  be  expected  since  tlie  t'./o  profiles  are  similar. 
By  comparing  with  cicample  1,  the  main  difference  in  this  ./ind 

T 

profile  is  that  the  maiiirnuiu  values  of  c.   occur  at  —  =  1.5   for 
3>  -^  q  =  o,   ami  also  the  band  of  unstable  wavelengths  shifts  to 
slightly  shorter  -avelengths  than  for  the  previous  './^sic  current. 
Moreover,  ".;hcn  u   =   -.000  km,  it  is  seen  that  ali   a/es  are  stable 
in  the  last  tv/o  cases.   Hencqi  it  may  be  concluded  that  the  stabil- 
iz;ing  effects  of  £     and  q   in  this  basic  current  are  much  more 
pronouaceci  v.'lth  increasing  D   than  in  example  1 
E::ample  3.   Doubie-jet  !  asic  curieiiL. 

In  preceding  t\;0  examples,  a  stai.ility  study  ■..■as  given  for 
syiTii^aetric  single~j_t  .■litid   profiles.   Now'  in  this  example,  .;e  ./lii 
investigate  the  stai>ility  characteristics  of  the  symmetric  double- 
jet  \/ind  profile,  represented  l^y  the  expressioii 

u=-6  (5in^v  -H^iH^y)  (2.16) 

'..-■hich  is  sho\m  as  citrve  (2)  in  figure  2. 

Inserting  (2.1C)  iiito  (2.12)  anu  performiUi^  the  series  of 
computations  yields  the  results  shown  in  figures  11  and  12,  cor- 
responding to  D  =  2000  and  4000  km,  respectively. 

Contrary  to  the  results  of  example  1,  it  is  seen  that  £ 
and  q  do  not  sIio;/  the  stabilizing  influence  for  all  wavelengths. 
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For  relatively  longer  wavelengths  (longer  than  the  wavelength  at 
uhich  the  iriaxlmU;;!  of  c^.   for  ^'  -  q  =  o   occurs),  the  d-af.f.ect 
iias  a  uestabili,::;.Lu}i  influence  ou  the  perturbations  in  Socli  a  ..ay 
that  the  magnitude  of  c.   is^  greater  for  19^0     than  for  B  =  o. 
Also  the  i.anci  of  v/aveiengths  for  which  instability  occurs  is 
broader  for  ^S  ^  Q    -      Furthermore,  for  these  longer  wavelengths 
the  effect  of  q,  a  measure  of  divergence,  is  to  further  increase 
the  instability  of  the  x/aves  similarly  to  the  ^—effect. 

It  should  be  noted  that  tv/o  maxima  of  c.   occurs  for 
D  =  4000  kra,   whe  .  .-S  ^  0  a^ic  9  4^  0   ;  and  that  the  v;aves  are 
much  raore  unsta'^le  v-hen  the  meridional  w'icth  u      1.-.-    i;icreaseil. 
Hence  th.e  divergent,  one-parameter  model  in  the  double- jet  v/ind 
profile  api-ears  more  unsta:  le  th.an  the  non-divergent  model,  v.'hich 
is  opposiLu  Lo    u..ie  results  for  the  single-jec  pro.Jiie. 

It  should  be  mentioned  here  that  the  subdivision  of  the 
r.asic  channel  by  finite  differences  into   6  mesh  lengths  cor- 
responding to  the  case  Iv  =  5,  is  prooably  a  cruder  representation 
for  the  dou:)le-jet  basic  current  than  for  the  single- jet;  hence 
the  results  for  the  former  may  he   expected  to  he    less  accLirate 
dL;e  to  larger  truncation  errors.   This  aspect  •..■ill  be  considered 
a:j,ain  in  i-'art  II. 

Example  4.     (j  _  g  ^  ^^.  ^  j^  _  ^^^  ^  ^;0  ^  (2.17) 

This  profile  has  easterly  v;ind  components  to  the  north  and 
south  of  a  belt  of  'i.^esterlies  in  the  center  of  the  channel, 
(figure  33).   Tb.e  values  of  c.   as  a  function  of  the  zonal  '..'a\-e- 
length  in  the  series  of  computations  are  sho./n  in  figures  13, 
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1'4,  and  i;^,  corresponding  to  D  =  2000,  4000  aa^-  6000  km^respecLivcly 

It  is  seen  tbat  ^  Uiux  4  \\d,-'^  a  j;.a  j.ii,;ing  influence  ay 
in  example  1,  but  au  uasta'.. le  Laud  of  ■>/aveleii;;jths  does  exist  for 
the  case   -^rzr  j^^jj.  and  q^=  1.3  >;"  10  ~  m  ~,  even  for  D  =  6000  km. 


The  rnaxiiiiuin  values  of  c.  (lioi;  B  =  <\   -   o)    occur  nearly  at 
L  =  1.0  and  are  greater  than  for  any  of  the  profiles  previously 

Li 

discossed.   However,  the  land  of  ..ave lengths  in  the  case  of 

^  ~  q  =  o,  for  which  instal^ility  occurs,  is  much  aarrov/er  than 
the  others . 

It  may  ;,e  concluded  tliat  this  profile  in  the  divergent  model 
:I.s  more  unstable  than  the  single- jet  profile  in  e::ample  1. 


^  u 


Fig.  16.   The  r:;onal  -...ind  channel  shov/ing  the  ci^ht 
sul:divisions  . 

The  '.eneral  validity  of  the  conclusions  ./hich  have  been 

readied  in  preceding  parts  may  be  q;.iestioned  recause  of  the 

finite-difference  niesh-sii'.e  corresponding  to  N  =  5. 

As  a  check  on  the  previous  results,  it  v;ill  ;:e  of  interest  to 

study  the  case,  K  =  7,  as  depicted  in  figure  IG.   In  this  case 

the  frequency  eq>-ation  (2.12)  is  replaced  by  a  fourth-degree 
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equation  of  the  form 

The  notations  in   (2.18)  have  the  foilov/ing  deflnitiouo: 
^/c^  K  -U,;  r  ,(//  =  .:,/  ,/i,  3.") 

•^  =.  r  +  d'f ,  . 

So=f:(T:'l) 

\4  zr-TcU-hD,  )-T,   , 

Vi  —   L'c/  ^'  I  I  ?     . 


(2.18) 
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Solutions  oT  eq.atLo.i  (2.18)  liave  been  oLLaiued  by  the  Muller 
metl-.od  [  &]  for  tho  same  series  of  examples  as  for  the  case,  N  =  j. 

Since  there  exist  differences  in  tlie  instai^ility  char- 
acteristics of  single- jet  and  double- jet  prof ilea,  it  is  of 
special  interest  to  examine  these  cases  using  the  finer  mesh- 
size  corresponding  to  N  =  7. 

Example  1,   Siiigle-jet  basic  current  (fig.  2,    curve  (1)  ), 
U"  =  h  (I  -COS^V)  (Same  a;^  (2.13)   ). 

The  values  of  c.   as  a  function  of  .vaveientli  in  the  various 
scries  of  computations  are  shown  in  figures  17,  IS,  and  19, 
corresponding  to  D  =  2000,  4000,  and  5000  km,  respectively.   A 
comparison  to  figures  3,  3,  and  6  shows  the  t".;o  cases,  N  =  3 
and  N  =  7,  to  i:e  very  similar.   Ho'.;ever,  some  differences  do 
exist  for  N  --  7  as  follo\/s: 

a)  iS   has  a  destabilising  influence  for  relatively  shorter 
v/ave  lengths  ; 

b)  tht  band  of  uiistahle  wavelengths  is  slightly  uider; 

c)  the  maximum  values  of  c.   for  jf?  =  q  =  o  occurs  at 

L  i  2.5  instead  of  at  2.0; 
D 

d)  due  to  the  lesser  stabilizing  effects  of   ^  aiid  q  \/ith 

increasing  meridional  width,   D,  the  unstal.ile  curve  (5) 
In  figure  18  for  D  =  4000  km  ano  curves  (2)  and  (3) 
for  D  =  6000  km  liave  appeared  (figure  19). 
Example  2.   Double- jet  basic  current  (fig.  2,  crve  (2)   ), 

U=B(Sin?y-HSir^^^)    (Sa..e  as  (2.1G)   ). 
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A3  indicated  earlier  greater  differences  miglii:  '.veil  have   co;i 
expected  in  this  case  because  the  iireater  coiiiiJlc;iity  of  the  dou  ie- 
jet  profile  is  not  as  \jell  represented  by  the  subdivision  into   6 
parts  (N  =  5) . 

The  resLiits  for  M  =  7  in  this  profile,  ii,i,u.:u   m  fi.-ures 
20  and  21,  are  also  quite  siuiilar  to  those  for  II  =  5  (fi,;ure3 
11  and  12)  except  as  follo.^s: 

1)  when  D  -^   2000  hn^; 

a)  the   aiaguitude   of  c.      is   nearly  double; 

i. 

b)  since  the  destabilising  influences  of  3     and  q  are 
much  greater,  th-j  aau  of  unstable  ..a-.  eie.i^Lns  nas 
v.'iuened ; 

2)  wnen   D  =  4000  km; 

a)  t  le   raar-ciiP.Uui  value   of   c.       for     fi    =  q   =   o   occurs    at 

L      ^      1.-; 
b 

b)  tiie  last  three  curves  (3,  4  and  S)  of  figure  21  no\/  sho  ; 

a  stabilii;ing  effect  ^/ith  increas.ln.;  q   for  v/avelen^ths 
given  in  this  fi;-iure,  but  it  seems  probable  tnat  the 
effect  of  q  '.;ill  le  reversed  for  sufficiently  longer 
v.'avelen;- ths  ; 

c)  for   P  =  fi^^  ana  q   =  l.li  x  10    la   ,   c.   vanishes 
at  L  =  :>000  hivi. 

In  summary,  it  may  be  stated  that  the  essential  features 

of  the  stabilit>  characteristics  are  similar  for  tlie  cases  N  =  5 

anu  7  ;   however  some  details  appear  to  \ c    lost  in  the  cruder 

subdivision  of  the  1  asic  channel  corresponding^  to  the  case,  N  =  5. 
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C.    Fourier-seriGo  mcthoc.l. 

If  a  series  of  tri^orc netric  terms  is  used  Lo  represent 
U  (y)   it  appears  appropriate  to  seeK  a  solution  of  (2.&)  for 
P  (y)   as  a  finite  Fourier  series  of  the  form 

4^(^3-21  ^y.  Sin  (71  >vV)  ,     n=/,  3  ,  .^ M         (2.1:) 

'■-'aere      X  — 

V      ' 

Then  the  expression  (2,15)  ouviousiy  satisfies  the  houndary  condi- 
tions at  the  ri[;id  '.7alLs  at  y  =  o  and  at  y   =-.  jj . 

As  the  valce  of  U   in  (2.19)  increases,  more  reliable  results 
regarding  the  staoility  features  may  !>e  expected.   V/iin-Nielsen  [3] 
has  found  the  sta'^ility  characteristics  for  -/a lues  of  n  up  to   3. 

Larger  values  of  N  (up  to  n  :=  8)   have  been  considered  by 
Eliasen  [Z]    but  mostly  in  cases  where  t'ne  effects  of  S      and 
divergence  terms  'vjere  excluded  by  assumptiou.   H'.s  results  sho„ec 
the  existence  of  only  one  unstable  mocie  corresponding  to  odd 
valines  of  n,  an  J  also  that  c.   at  ~  =  2.0   is  the  same  for  3-Yl-^, 
but  some. 'hat  smaller  for   a  =  7  aiid  8.   Hence  this  study  v/ill  oe 
concerned  only  .;ith  small  and  odd  values  of  n. 

To  verify  further  the  results  in  tbic  foregoi.ig  section,  the 
symmetric,  hariuonic  ;;onal  velocity  profiles  uiil  again  be  considered. 

1).   Single- jet  basic  current: 

U  =  B(/-CoS^Y)         (same  as  (2.13)   )     (2.20) 
'..'here  B  =  30  m  sec. 

Considering  the  case  W  =  5,  inserting  (2.20)  and  (2.13)  into 
(2.6)  yields  a  system  of  linear  homogeneous  equations 
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[|(^-3X^)]t+-[(^^'-^i'^+T;;Cfp-B(^>^^^  (2.2i: 

lu  order   to   have   uon-trivial   solutions   for   0    ,      ^.  ,      aad  0^    , 
tiie  determinauL  must  vanish,    ./hicb  ,-;ives    tlie   frequency  equation 

C'CRWE)-  clEvR\/+SUi)^RLii^]+c[bf;^V-TH)f^^R^'iSu3) 


-PHivn^[&(SV-TH)-PHVJ     =-0 

T  --■-  |-  ( 2  /  A"+  // '  )  ,       U'  =  Z^  A'-^-P "  -h  Q ^       ; 

and        V    -=^  p~d  (^^A-+M^)     , 

In  the  case,  l-l   ~   7,  t~ne  frequency  equation  (2.22)  is  replaceu 
by  the  fourth-degree  equation: 

C'^CYRWt  )tC'(YK.f.KKWt:')-f-C'(YK,fZK,-TRFx;-hC[YK^ 
+2K, -(ES+-Ke)TXl  +[2K2-((rTS-FH.  )TXJ  =-0 

• /her c  X  =  -f-  (  /./"  4- ^:^  K'  )  ,    Y  -  ^9  X  +''-^^'     , 

M.^SV-TH  .     Ko-Eric  +  RUirT     , 

K,  =Eri, +CTno-?Hu.s 

and         K,:3  6-n,-FHV         . 


(2.22) 


(2.23) 
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The  roots  ol;  (2.22)  ant!  (2.23)  for  the  previous  series  of 
computations  are  sho./n  in  figures  22,  23,  24,  25,  26,  aa...  2  7, 
corresponding  to  D  =  2000,  4000,  and  6000  I;tn  for  tl.e  cases, 
K  =  J ,   aiici  N  -   7 . 

There  is  O'viously  a  close  similarity  in  the  essential 
features  of  tlie  stability  characteristics  anoa^  the  cases,  N  =  3 
(  jivea  in  ki/iin-iUeLsen  [-.j)  ,   i-j  =  5  anc'   N  -  7.   Fur then.iore, 
tliere  also  exist  ciuillarities  i  etu-e-^.i  '--i'-.  .■u^l.j.l.c.  i...  l..l.  _^.ij.tc^- 
difference  iMethoci  and  those  of  the  Fourier-series  metho...    'X 
may  he  noted  that  the  results  of  the  case,   d  =  —  (M  =  5),  in 
the  former  method  are  particularly  close  to  luc  r'_o..its  of  th.e 
case  N  =  3   in  the  Fourier-series  raetlioc  ,  ./hile  the  results  of 
thiu  case,   a   =  —  (N  =  7),  correspona  more  nearii  to  those 
for   h  =  J. 

It  may  !  e  concluded  that  the  results  for  sinall  values  of  N 
jenerally  represent  th.e  essential  features  of- the  stability 
cliaracteristics  in  tlic  quasi-i.arotropic  flov;;  ho'./ever  the  foilo.;- 
ing  differences  '/ith  increasing  N  do  appear: 

a)  the  rnaxinrain  value  of  c.   for  p?  =  q   =:  o  is  less 

(c.     =  9.57  rn  sec   for  '^   -   j   ar^o  i3.40  m  sec 
1  max. 

for  N  =  7  ); 

b)  the   sta''ili;.in;i   effects   of  &  and     q     v/ith   increasing     D 
are   less   proiiounceu; 

c)  the  niaximuvii  \alues  of  c.  are  shifted  to'.'ard  longer 
v/avelen^ith.s  sich  that  for  the  case  ^  =  o  ano  q  =  o 
they   occur  at  -  -  2.0,    for  W   =  3;   -  =  2,5,    for  N   =  5; 
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and  —  •-  3.0,  for  K  =    I; 
ri)   thf  '.-and  of  iional  \/ave lengths  for  ■'./hich  instal  iLity 

occurs  becoines  ..'ider. 
2).   Douijle-jet  basic  current: 

U  =  B  C/-  CoS^y)  (2.V4) 

where  B  =  30  m  sec 
This  profile  has  U  =  o  at  both  rigid  /alls  and  at  the  center 
of  the  channel,  a.id  also  lias  U  =  2B  at  y  =  r  anv  at  y  =  7  D 
(fig.  34,  U   (o)   ). 

By  inserting  (2.24)  aviv  v.:.19)  (K  =  5)  into  •,-..),  the  foilo- 
ing  set  of  linear  homogeneous  equations  is  o'  ta L  .   : 

[f  (/^A"-//-)]  i,  ^-[(•9/W/''"fq^)C-f-p-p(9-v4-y^'  )  ]  4».  =0    (2.23) 

In  oruer  to  have  non-:.;ero  solutions  for  ?!- ,  /   aa;'  if    ,    the 
deterniinan'c  must  ./anish.   This  couu  i  Liur;  <j,j.ves  the  frequency 
equation 

where  t  -^X+y^'^-i-U"       ,         ix  =^  p  ~  B  (^  "^  i^' )     , 

H  :.:  E(^(>c-un         ,      T  =  4  (9A:'  -f //"  )  , 


(2.26) 


^ 


=^  8  -H{^Ar+-jr)  ,    {^j^^^^f-^  i^'^T, 


:i     \/  ^  :3  -p  (^^x:-f-,//-)   , 


anc   ,,  _  ,,  _.  ^  ^^ 


This  frequency  equation  (2.26)  has  been  solved  numerically 
for  the  previous  series  of  coinijutatioiis .   The  values  of  c.   ./ith 
respect  to  ivaveleugth  are  shovm  in  figures  28,  V)   aud  30. 

It  is  seen  tliat   ^  aacl  q   effects  '..o  not  give  a  stai/i lining 
influence  for  all  wavelengths.   Fur  relatively  longer  v/ave- 
iengths,  ^     and   q  liave  a  destabili;iing  influence  on  the  pertur- 
bations ir.  the  sense  that  the  band  of  unstar.)ie  './avelengths  is 
much  broacier  tnan  in  the  case  R   -   o      or  q  =  o  .   Moreover,  the 
.7aves  are  nore  unstable  v/ith  increasing  meridional  '-/idth  D. 
Hence  the  divergent,  one-parameter  model  in  thi::  'douiile- jet  v/ind 
profile  is  mere  unstable  than  the  non-divergent  model.   Tliese 
results  are  consistent  v;ith  those  in  th.e  previous  example  for  a 
double-jet  profile  in  '.fnich  the  method  of  fiaite  differences  has 
been  used. 

For  tliis  profile,  if  even  values  of  u  (n  =2,  4,  6  )  in 
(2.19)  are  included,  there  exists  other  amplifying  raodes  different 
from  those  previously  discussea,  v/hich  are  not  shovn  in  this  paper 
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3.    Stability  properties  of  the  stratified  moJel. 

The  vorticity  equation  for  this  model  has  ;>e(jvi  outaiaed  by 
'Araason  [8]  and  may  l:.e  v/ritten  in  the  form 


( V  ^-  T ' )  If  +  V '  V '  V '  ¥^ + n  -  g '  V  •  vv^  -  0 


(3.1) 


where  q  ^^  ^jr^^     ^    V  -  9<  X  V  f  ,  ^^'^   V  =  U  i  . 
Here   H   is  the  thickness  of  the  stratified  layer,  aad  j.  '   is  the 
acceleration  of  gravity  corrected  for  compressi  iiity. 

By  using  the  pertun^ation  method,  (3.1)  ca   e  reduced  to  the 
form 

where  "^    is  the  -perturbation  stream  LL^action  and  U  =  U  (y)  is 
the  mean  zonal  velocity. 

substituti:-    Y'(x,\/.  /■)  ^  c|>c^).e^^^^^"^'^^ 

into  (3.2)  yields 

^•^-c)(^"-(r<t')+-^s-DVq'(c--D  )]  |.  ^o  ^.s) 


where  the  primes  c;eaote  differentiation  x/ith  respect  to   y. 

In  comparison  with  (2.6),  equation  (3.3)  has  one  more  term 
v;hich  is  peculiar  to  this  model.   Thus  it  may  e  exfjected  that  the 
stability  characteristics  in  this  model  "./ill  differ  somev/hat  from 
those  of  tlie  previous  divergent,  one-parameter  model. 

Th.e  analysis  v.'ill  be  restricted  to  the  initial  ;:onal  iJind 
profile  of  the  type 

U  -  B  C/-  C0.^^  ^)  (same  as   (2.13)   ), 
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aiici  also  to  the   casic  N  -  J   in  th^  expression  (2.15).   The  result- 
ing frequency  eq^at i.r-.  : 

C^(Ru)E)  +  ^--[E(Rv^SU))  -f-RUJ(r]  +C[E(SV-TH) 
+e(RV+^W)-?HuO-hi&(SV-TH)-fHV]  =  0       ^"''^ 

where  E  -  X  +//•'  +  q'     ,   G  =  ^  -  f- {  3jf^3f-K  )  , 

and  v=  ^-B  (2i\^  +  ir-f-  qM 

The  results  oi:  the  previous  computation  series  are  given  in  f injures 
31  and  32,  corresponding  to  u  =  2000  and  4000  km. 

It  is  seen  by  comparing  the  curves  (2)  and  (3)  iu  fi^^ure  31 
that   q   has  a  aestaiviliziiig  influence  Cor  reiati-.ely  longer 
'..'a -e lengths  (loii.^or  than   the  waveleagtn  at  -/hich  the  maximum  c, 
for  curve  (2)  occurs).   Moreover,  these  effects  are  verified  v.'itli 
increasing  values  of  q,   as  given  by  curves  (4)  aiid  (i3).   Hence 
the  divergent,  stratified  model  is  more  unstable  thian  the  non- 
divergent  one.   Also  comparison  of  figures  31  anu  22  shous  that 
the  stratified  model  is  more  unstable  than  the  previou:^  one- 
parameter  model. 

These  results  are  also  verified  by  figure  32,  \.hich  contains 
similar  results  for  D  =  4000  km.   Note  tliat  t'ne  .;aves  are  more 

stable  with  increasing  D.   Nevertheless,  for  the  case  ^  =^  .-^/fj-   , 

o  -12-2 

q"  =  1.3  X  10   m    and  i)   =  4000  km,  it  seems  possible  to  have 
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a  I  and  of  unstable  './aves  at  loager  .;avelen^ths  which  are  not  sho'.m 
in  this  figure. 
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4.   Momeatuiii  changes  •■ 

It  is  of  special  interest  to  see  ..Viat  i.iay  .  e  iaferred  from 
the  solutions  of  this  particular  problem  about  the  effects  of  the 
unstable  perturl.atious  on  the  basic  flov/  through  momentuiu  transfer. 
Changes  iii  the  :ional  current  './ould  alter  the  profile  so  raciically 
that  the  initial  solution  of  the  perturbation  problem  no  longer 
applies.   This  implies  that  the  basic  profile  may,  in  turn,  alter 
the  structure  oi.    tuu  perturbations. 

Recently  i.i  the  paper  by  Wiin-Nielsen  [3]  ,  it  has  i:een  found 
fnat  a  single- jet  zonal  v;ind  profile  initially  can  change  to  a 
LOU!>i.e-jec  prociie  in  one  or  t./o  uays  aat;  to  meridional  momeutuiu 
transport  by  tl^e  eddies. 

It  is  obviously  of  interest  to  determine  hov;  an  initial  double- 
jet  wind  profile  ':/ill  change  as  a  res-..lt  of  momentuia  transport  '.■y 
the  unsta'ile  pertur  at  ions. 

In  tb.is  section,  we  will  consider  the  double- jet  i)asic 
current 

JJ^B  U  -CO^^^y)  (same  as   (2.24)   ), 

■..'hich  is  designateu  iy  U  (o)  in  figure  34. 

From  the  equatioii  of  motion,  the  changes  in  the  mean  zonal  current 

can  be  expressed  in  the  form 


.   ^.=  -  ^'  •   (4.1) 

at       a  ^ 

where  the  prime  designates  the  perturbations  and  the  bar  refers 
to  the  zonal  mean 

u  ^ 

26 


Thus    the   cha.'..:^;,.  .id   profile   can   be 

compuLc  ^L'.ji.i  L.i;c  i..^i.iui(..'i,ax  co.iv'wj-'^ti.icc  of  the  ->- '-y  .  lu.i'.entura 
transport . 

Ii:    order    to   evaluate    this   quantity   fron  the   solutio.i   cjf    the 

jjert.^r    ai-.^u;:    ^.r:-)..'i.ijrn,    ticiuaun.  :    l^-.i;    :viay    .e    ..ritten    ia    the    form 

Tlie  notations  have   the   folLo-'ing   definitions: 

4>,'  ^  Oi  ^  I  o.i 

tr=  /•-'  +-i  '^^~ 
...lere    tiie   superscripts   denote    the   auiplitudes   of   the   perturuat^oas 

(see  equation  {l,lj)      )   corresponding  to  the  t'.;c  roots  of  the 
frequency  equation  (a  quadratic  eqi;.ation)  ,  c    a.^-J  c^  , 
respectively . 

A  derivation  o:^  this  equation  has  been  given  by  Wiin-Mielsen 
[  ,   .,^  .  -v^..  ~1  .   llcirevor,  equations  (3.C)  a.;-.-  (3.10)  in  his 
paper  must  be  rey laced  by 


ar.'j 


G=f(7^--/)      .  (4.5) 
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—  /    ■ ' , 


(4.6) 


conditioiis,    (4,2)  may   'oc  v/ritten   in   tae?    form 


inte;;rat-(.ori, 


—  I 


(4.7) 


1  -f   .  O  y 


.(^-:) 


:cj.l-ii.--      u 


rariit 


froii.   (4.^)'vhat    jIII. 


;:   .:iay  ea 


seen 


■  lal   r'.onal    ..'in' 


i>ro£iie   wich   ti.ae.      ilia   iunc;:^ 
-■iiitive  maxiinuia  at   tlie-  ceater   - 


(cosh.  (2  //  c  .  t)  "1  ,' 


i  i  e  ir  c3 1:  o  r  ' 


proi;i  j.e 


■J(o) 


;y),    ijiiO./u    iu    :.:i,5ure    Ji,    has    a 
:L,    au-J     .•  ,   ' ::  ^  - 

Fuvtbex'mor e ,    ;■■  inc c 

^o.    I''-    ■  .  '..ras    thau   tiie 

as    the   sc  ) . 

',iii,  1:101  Kriiiimurn   O-    '.        .ci.al 
t.    Lue   center   of    the   cha.ti.ei   v;ili    Leud    Lo    increase 
.    .iiitually  the  maxinran;  .-         .-O  appear  at 


:       0   . 

■3i: 

^  -  V  S  Jl  C  J..V  t; 

It   ai^ 

:.    that 

t  ' 

tijG   double-jeu  uivid   jjrp.l'iie   will   tenci    to   cliaa-e 


1  J.  e  . 


Tliis  aieans    i,n; 
i:iLO  a  siiv;;ls 

:t    value's    :;-:      D      a.: 
tioa  )      as   a   function   of 


urc    V- 


.••4   ooi. iiai^iii 


;)r  noi: 


;  ■-   .3000  i;;-i  aiu!  a  .,oaal  ..avele;i:; th 
1,   =  2000  km.   The  .-a-.e  iyuasc  velocity   c    ..    Laiuec;  for  'i:»:e 
case     -  3   in  (2.1Q;  e  fi-ure.  The  ./ave- 

c-^t'^  of  2000  I;m  cojTi -St^^w.-ed  to  the  maximum  iastaM'ii  y  '-.'hen 
D  ^  3000  km  (in  this  case   c  =  28,9  ^   i  3.9  m  sec  '" ) . 
Ti:e  am'^;litu':'.e  k-.     :>'^    the   ";ave  component  ovie  of  the  ix^itial  pertur- 

acj.oii  ^Lj.ua.a  .-.■CLj.^j-i    .O.C.    i,ii:'c.^-.i   to  correspond  to  a  meridional 

ind  of   5  m  sec    (A^=  1.59  x  lo''  m  sec"  ),   hile  A^  =  0,  auc 
E  =  30  m  sec"*",   lu  general,  A^  and  A,,   are  related  to  ;i'y  as 

.oiio-./s: 

r  -  h     =-  ^^    . 

After  2  1/4  days  the  initial  profile  has  chan^^ed  to  the  siagle- 

Figure  34  shoulc;  ;..e  compared  with  figure  ^.   :.ere  the  para- 
icters  iiave  the  same  values  except  that  A^  ==  A,^  =  1.59  x  10  m 
:c'^'.   A-ain  a  :^f:r.,^  ,,  of  the  double- jet  prof ile  iato  the  sin-le- 
'et  profile  has  occurred,  hut  at  a  sli;^htly  faster  rate  for 
...,  \   o  tiian  in  the  case  A,  =  o. 

i  .  ,ui.  ^  -.■..  contains  the  distruution  of  U(y,t)  as  a  function 
time  for  D  =  6000  km,  L  =  '^000  ..m  (corre^ponein^  to  maximur.. 
instability  for       'j)  ,    ano  other  parameters  as  in  figure  35. 
I.,  uhis  case  th:.  ::^^r^lw.,  rate  is  muc>  .lo;er  with  the  lar-er   U 
a^.d   L.   This  al"o  can  le  shown  by  examination  of  equation  (-i.'J). 

It  is  obvious  from  (4.-;)  that  if  the  imaginary  part  of  the 
phase  velocity  vanishes  (c.=  o),  .hich  means  Varotropic  stability, 


the  zonal  ■.-/inc  ^        is  independeat  of  time.   Hence  the  changes 

preuictea  in  the  r^onal  profile  £co:a   these  coiviputatio.is  are  i'uc 
to  the  presence  of  uiis table  barotropic  ./a.cG. 

Thus  it  may  I'e  coiicludecl  that  if  there  exists  'arotropi-c 
in;-; tar. ill ty  in  a  lasic  flo',-  cliaractcri^.e!.'  'y  a  ..•oul  ie  maximum, 
it  is  a  dynamical  coasequeace  of  the  couvergeuce  of  meridionai 
transport  of  moinentuva  hy  the  eddies  that  the  ;^oii.al  current  tends 
to  develop  a  single  maximum  of  zonal  ''/ind  spcjCu  .   This  results 
from  momentum  being  concentrated  in  the  center  of  tlie  channel  by 
the  eddies  at  the  expense  of  both  maxima. 
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'j .        Summary. 

In  secLica  z  tiie  staL-iiity  of.'  ;.aro tropic  lio.:  lior  the  diver- 
gent, one-pai*ametei-  mociel.  has  been  investigated  for  various 
syirimetrical  :ioual  v^in-^.  fields,  including  sin^^le-jet  and  double- 
jCt  profile.s.   'i''.;o  uiiferent  ap^jroxiroatioas  arc  used  for  the  ampli- 
tude of  the  perturbation  stream  functiou,  one  of  chich  is  based 
ou  a  method  of  finite  differences,  and  the  seconc  on  a  represeuta- 
tion  iy  a  finite  Fourier  series. 

The  i.asic  equation  in  this  study  is  the  vorticity  equation 
'.'hichj  for  a  nar-'.-o'ilc  v'ave  solutinn.,  reducer.  i:n  a.   secoufi-orcier 
differential  equation  vita  variable,  coefficieats  for  tl^e  perturlja- 
tion  ampiitude. 

The  main  coiiclus  r  .•  •  -  :^re  that  n-'  ■■    'o  i:ie;:hod'-  ;^ive  essen- 
tially the  saiiK  .-tcoudly-j  uue  iutroduction  of  divergence 
into  the  model  in   the  form  of  a  Helmhcit<i  term  tends  to  reuuce 
the  iubta-  ility  presc      l.  .,  ..j,;- divergent  .u^  .^^  for  a  single- 
jet  zonal  wind  profile,  '..tiiie  the  reverse  holds  for  a  double- jet 
profile.  Thirdly,  it  is  also  found  that  the  \/ave3  in  the  ciiver- 
lier.t  model  are  rriorc  stabu,  ..lui;  i,  .creasing  ineriai.^-ui-^  ./idth 
betv;een  the  assumed  ri^^id  boundaries  for  a  .;i;.jle-jet  profile, 
'./hile  the  reverse  a^ain  Isolds  for  a  'oul;le     profile. 

Secti'ju  J  is  concerned  /ith  tbe  sl        j.ia^aci...-.  i;,  lIcs 
of  'Aruasou's  "stratified"  Uarotropic  model  ..     a  finite 
Fourier  series  for  the  perturi-ation  stream  functio..   If   a-:- 
foL.ad  t.iat  i-'he  '..  ivtrv^fni. ,  stratified  model  is  more  unsiiai-ie  Lhaii 
the  non-divergent  pue;  and  also  even  more  unstable  than  the 
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divergent,  oae-pax-amtjCer  mociei. 

In  sectioii  -i  ^..,^   second-oroer  effecti:  w.  ;.'ae  unstauLe  perUur 
datious  on  an  initial  double- jet  zonal  fio;/  are  iavestigated , 
It  >'as  JJouHi'  tliat  this  type  of  profile  of  the  :'.0:ial  wind  ch^a^es 
to  a  siii'jie- j^  ^   :j..al  ^.^o..  '.;.  ..    .  •  a^  ■-. ,  solely  ...^..e  Lo  the 
convergence  of  meridional  transport  of  niomentum  in  tlie  center  of 
tl).c  channel  caused  by  the  unstable  barotropic  perturbations. 
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Fig.  4.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.12) 
as  a  function  of  v/aveleagth  for  the  same  parameters  as  ia  fig.  3 
except  D  =  3000  km.   The  curves  as  in  fig.  3. 


35 


Tj                  1   _M.  -  SlL.S   _     '     _U.1_. 

1   1   1   i       1       ;                    1                                        ;    1    1    1 

III 

'              ■                                                            '    '    '    : 

III 

^              '_,;"""..] +    h  -t    ' 

!  :     i     1                                           !  :  :  i 

:  '      1 

si                          i-    ■  -^                         I^ 

'  i  '     1     1                               1  i          i  ■  '  ' 

1      1 

?^                                                                           M                                                                  ! 

1 

v^                                    ^                i    ! 

i    i       :        i        1                                                    1                   •    : 

1 

1 

.^                                             i        •_ i         I 

i  j        '    ■■  '      1                               1    : 

^  i  i  1 

*♦■'                  1                                   ill 

'   i  ^'    '  ' .                                1                   11; 

'  1 

1    1    '            i                            i 

\>^        ,  i  ^<.,                                 1  1  ; 

:  i  1 

H            1                              III                     -r  . 

>''^i  '■     1  1     .         "v.  ^'     ■     1           '  !  :  ! 

■  . 

1                        1/ 

1       1    I       1       1              i       1      ^/  J        1       1    i           1    I    ! 

'  1     ' 

/ 

1  1 

1 

/ 

>^. 

i  1 

/ 

\      -r-       -L 

1 

/ 

X       I                1 

1  i 

/ 

1     X     1  1  i  1 

f 

1          jS    1   1   1   1 

it     lt7 

^    1    1    1 

1 

T 

IN    !    ' 

1  1 

6                                       / 

i  \ 

!   1      1 

r    ± 

i  1  i\ 

1   '      1 

T  _ 

j  i  '  \ 

■   1      i 

/'    >  — 

1     <     1 

\.i   1   1   1 

^  v/   ^ 

1    j    ' 

N.'    1    ' 

/  /  '             \j 

\  1 

/  /                  \ 

1                                                          i  i  :  , 

,     \ 

/  /  !                   > 

.,  j      ,        ,                                             \  \  I 

:    \i 

/  /                        ^ 

i                                           1                 ill' 

i    '    ;  \ 

11 

\C2)                                                          '  '  ' 

i    1    1^ 

4                                                   '/                         1 

\   ^\           '                             1                   i  1  ' 

•II' 

1                          j  /A 

i\'  1  1                           '                                ill 

. 

1 //  ] 

\  '  '                                                           III' 

1 

"7/  f"^-^^ 

;   i   1   i 

//  /'^        \      ^~ 

T~  ~^                             "   n  {~i~' 

1 

-r 

1 1  \/            if  ■> 

\ 

'   I 

1 

//'  /                    ^^v' 

\                                              ' 

I   1 

III                            \ 

\ '                                            1 

'   1 

1 

//  /                           \ 

\i 

1 

2                                          //  /                          \ 

\                                                                                    1        1        :    1 

1   1   1 

//  /       '                '   '  \ 

\        1                                                   1                                  •        ,     ' 

;    i    ; 

/  / /     /^  1                   \ 

i                                                                                                    '    ' 

i  1 

1 

1 //  /     A  1  Y                   \ 

\                                                                                                '    ' 

j 

1 

7  /  /  /'    /\y  Y          \ 

1 

ill/        ^  v^             '  \ 

1  1                                                                                     i    i    ;    ' 

1 

III               \   '                  \ 

1                                                                                    I'll 

' 

1 

1 1 1  1           v_             \ 

1!                                                                                                                      i      1      '      : 

i    , 

1 

///  /         iV          11 

Ij                                                                                                                      1      '      ' 

'  i  I 

/  /  /   /  1            y              1  \ 

11                                                                                                                      I'll 

I'll 

0                               / /  /  i             \i              il 

ill!                           1                 i  '  1  : 

'        1    1 

10 


12  L(io^\<m') 


Fig.  5.   The  Imaginary  part,  c,  of  the  solution  to  eq  .ation  (2.12) 
as  a  function  of  vjavelength  for  the  same  parameters  as  in  fi^.  3 
except  D  =  4000  km.   The  curves  as  in  fig.  3. 
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Fig.  6.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.12) 
as  a  f  aiction  of  Vi/aveleagth  for  the  same  parameters  as  in  fig.  3 
except  E  =  6000  km.   The  curve  (1)  as  in  fig.  3. 
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Fig.  7.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.12) 
as  a  function  of  the  meridional  width  D   for  L  =  5000  km  and  the 
same  parameters  as  in  fig.  3.   The  curves  as  in  fig.  3. 
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Fig.  8.   The  ima^^inary  parL,  c .  ,  of  the  solution  to  equation 
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Fig.  9.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.12) 
as  a  function  of  wavelength  Tor  the  same  parameters  as  in  fig.  8 
except  D  =  4000  km.   The  curves  as  in  fig.  3. 
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Fig.  10.   The  imaginary  part,  c,  of  the  solution  to  eq  ation  (2.12) 
as  a  function  of  v/avelength  for  the  same  parameters  as  in  fig.  8 
except  D  =  6000  km.   The  curve  (1)  as  in  fig.  3. 
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Fig.    11.      The    imagir.ary   part,    c.,    of    trie    solution    to   eq.ation    (2.12) 
as   a   function  of  \7avelength    for     U  =  B(siu(  7ty/D)   +   sin(37ty/D)) , 


d  =  D/6,  D  =  2000  km,  and  B  =  40  m  sec" 


Tbe  curves  as  in  fit 


3. 


42 


10 


-2  L  (jo^krn) 


Fig.  12.   The  imaginary  part,  c.,  of  the  solution  to  equation  (2.12) 
as  a  function  of  wavelength  Tor  the  same  parameters  as  in  fig.  11 
except  D  =  4000  km.   The  curves  as  in  fig.  3. 
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except  D  =  4000  Km.   The  curves  as  in  fig.  3. 
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Fig.  15.   The  imaginary  part,  c,  of  the  solution  to  equatioii  (2.12) 
as  a  function  of  wavelength  for  the  same  parameters  as  in  fig.  13 
except   D  =  6000  km.   The  curves  as  in  fi^.  3, 
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Fig.    17.      The    Imaginary   part,    c,    of   the   solution   to   eqi.ation    (2.18) 
as   a    function   of  vvavelength   for     TJ  =  B(l-cos(2;Ty/D)) ,      d   =  D/8, 
D   =  2000  km,    and  B   =  30  m  sec"'   .      The   curves   as    in   fig.    3. 
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Fig.  13.   The  imaginary  part,  c,,  of  the  solution  to   equation  (2.18) 
as  a  function  of  wavelength  for  the  same  parameters  as  in  fig.  17 
except  D  =  4000  km.   The  curves  as  in  fig.  3. 
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Fig.  19.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.18) 
as  a  function  of  v/avelength  for  the  same  parameters  as  in  fig.  17 
except  D  =  6000  km.   The  curves  as  in  fig.  3. 
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Fig.  20.   The  imaginary  part,  c,  of  the  solution  to  equation  (2„18)  as  a  function 
of  v;avelength  for  IT  =  B(sin(7ry/D)  +  sin(37iy/D) ),   cl  =  D/8,  D  =  2000  km,  and 
B  =  40  m  sec"'  .   The  curves  as  in  fig,  3. 
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Fig.  21.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.18)  as  a  function 
of  wavelength  for  the  same  parameters  as  in  fig.  20  except   D  =  4000  km.   The 
curves  as  in  fig.  3. 
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Fig.  22.   The  imaginary  part,  c,  of  the  solution  to  equation  (2,22)  as  a  function 
of  wavelength  for  TJ  =  B(l-cos(  J7ry/D)) ,   N  =  5,   D  =  2000  km,  and  B  :=  30  m   sec"'  . 
The  curves  as  in  fig.  3. 
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Fig.  23.   The  imaginary  part,  c . ,  of  the  solution  to  equation  (2.22)  as  a  fuuctioa 
of  wavelength  for  the  same  parameters  as  in  fig.  22  except  D  =  4000  km.   The 
curves  as  in  fig.  3. 
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Fig.  24.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.22) 
as  a  fanction  of  v/ave length  for  the  same  parameters  as  in  fig.  22 
except   D  =  6000  km.   The  curves  (1)  and  (2)  as  in  fig.  3. 
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Fig.  25.   The  imaginary  part,  c,  of  the  solution  to  equation  (2,23)  as  a  function 
of  wavelength  for  the  same  parameters  as  in  fig.  22  except  N  =  7.   The  curves 
as  in  fig.  3. 
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Fig.  26.   The  imaginary  part,  c,  of  the  solution  to  equation  (2.23) 
as  a  function  of  wavelength  for  the  same  parameters  as  in  fig.  23 
except   M  =  7.   The  curves  as  in  fig.  3. 
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Fig.  27.   The  imaginary  part,  c,  of  tlie  solution  to  equation  (2.23) 
as  a  faiictioa  of  wavelength  for  the  same  parameters  as  in  fig.  24 
except  N  =  7.   The  curves  as  in  fig.  3. 
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Fig.  28.   The  imaginary  part,  c,  of  the  solution  to  equation  (2,26) 
as  a  function  of  v/avelength  for  IT  =  B(  1-cos  (4Ky/D) ) ,   N  =  5, 
D  =  2000  km,  and  B  =  30  m  sec"'  .   The  curves  as  in  fig.  3. 
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Fig.  29.   The  imaginary  part,  c,  of  the  solution  to  equation  (2,26)  as  a  fuiiction 
of  v/aveleugth  for  the  sarae  parameters  as  in  fig.  28  except  D  =  4000  km.   The 
curves  as  in  fig.  3. 


r^ •    1          1    1 

Mil      MM            !  1 

1         1         >         !  I   1 

1           i 

" 

'•0     ~"    it         T4I- 

M  !     i  M  ■            ' 

1                      M    ' 

i 

^               rx- 

M         i      M 

Mi             !             1    !    ^ 

i 

i                   \     zr 

1  i  1  1      MM 

M         ■    !     i  !  '  1 

1 

^                               !          t 

1       .1                    1  ■  i  '■    '  M  ' 

!    1           ]               ■    '    ;    ■ 

•       1 

;^        1   :                   i      ' 

1    ' 

—         1               i     i 

/                 ■  ^\'          1     !     1               i          M          '     i     1     ■ 

1           1           \      '■■■'•    ■ 

1    ' 

O         ' .      _^ !   /;  J, 

/       i  1  \  i  1  i          1      M  '  ' 

i         >         I'M 

1       1    1 

t/*"^        >y[      1                     M       1                                  1       !             j       '       1       ^ 

i        1        i      ,      M 

'    !  i 

10          "^               'Ty 

/  .      '  \  1      V  ■  M              M-      1  !  '  ' 

i        1  . 

'  1  1 

Zi^         ■\l    \'      '         '      •      '   '      '         i  '  ' 

i    1 

'                   // 

;     Mi 

1  1  1 

1 1  \ 

\     i  \                 1  i                       1      1 

1     1  '  1 

j    1 

1  L 

\      \                                               1 

1     'III 

1  1    1 

//    1 

V       \                                  '    '                 ' 

i    1 

//     ' 

\     \                        1      M 

111' 

.  1 

— yy    1 

\   i\                              ' 

1   1   1   : 

1    1 

__//    Jk 

1  \  j  \               1               1 

! 

7/     "^^ 

..—                         \  \              1 

i     i  ; 

i 

S                      //    X 

|f'"'\                     ^\             1 

i     • 

1  I 

'         \           1  1  '  1  >-i      M      i      ' 

1                       •            \            ■: 

1  1 

//         ! 

\|                                     1      \     I'^K^    1             j            M 

1                      111 

1    ! 

/!  /    / 

\,                     '     \      i^r~  ^      i    ' 

Mil 

— ip J/7X- 

-i-      -^          \    ^--^-T^--^^    - 

1    i 

- 

---]-L-L///--(^, 

;:-i— 1 — hTrM"4-t-ri — !    M    1    i  !  " 

- 

1     11 1    ^ 

Nv                  '     1        '\M  '      i      i  ■        i  '  1    . 

1    11 1  / 

i^v               i     -'IVJ^- — ' — 'I  \T 

1        !         1                 1                 II                               "*"««^ 

'■'    '    , 

^           li//     ' 

:    1    '    i               1    1        1    1 

i.-.^ 

~Hnf  I  ^~f 

-M-    ^^^  -"p^    ~Ht~^ — i  j  h    j  i  1  i- . 

— i--n^— .  .  ;■■-. —  1  ,  M 

1           1 

- 

//  "  f^ 

j                    Jjl---^  "^'"'^'^ 

i      '      !      1 

!        1 

i           i 

X   l| 

1  /  / 

s,  1                          -^-"^^     '    '  \ '            1 

i           1 

1  /  / 

v^      _t^»*  "*            M  I'       i 

!     1     1     ; 

1 

7  / 

1                           !     \i        1 

1 

1 1 

1                           1  i     ' 

1         1     ' 

[ 

1 1 

1                          ' 

MM 

1 

'^i-     1—                     /' 

1   1               1    ' 

1 

r 

1          \  1     1     M 

:             i 

1  1  1 

jl 

\  \     1     1  1 

i 

1  1  1 

|! 

'\!      !      ! 

!     1     1 

1 

i         1    1 

1 

'                '     i  1' 

i  ! 

1 

1 

1  Ij 

1  i  i 

i  V 

1    '    ! 

! 

!  V 

i  M  ' 

j              i 

^  t: 

;     1 

1           MM 

1         i     1 

9 

'                '                                '1 

1         1     i 

i                          111 

1                   1 

1         I     1 

i                          111 

i                                                 ! 

1        ! 

1          '1              i 

i 

1        1 

1        ;l      1        i 

1           : 

1        i 

1                                    1        '  I      ' 

1 

1        1 

1  1      ' 

i        j 

111 

'               i    ' 

'  !      i 

-1    i     I 

1        j        1 

1               i 

i  '      1 

- 

i        1        j    1 

1               I    ' 

1      ' 

0              1...^ 

1     !     !  '     : 

I'll 

!    1        1 

10 


12  ],(/o3  \crfi) 


Fig.  30.   The  imaginary  part,  c,  of  tb.e  solution  to  equation  (2.26) 
as  a  function  of  wavelength  for  the  same  parameters  as  in  fig.  28 
except  D  =  6000  km.   The  curves  as  in  fig.  3. 
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Fig.  31.  The  imaginary  part,  c,  of  the  solutioii  to  equation  (3.4)  as  a  furictioa 
of  wavelength  for  the  stratified  model,  TT  =  B(  1-cos  (27ty/D)),  II  =  5,  D  =  2000  km, 
and  B  =  30  m  sec"'  .   The  curves  as  in  fig.  3. 
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Fig.  32.   The  iuiasinary  part,  c,  of  the  solution  to  equation  (3,4) 
as  a  function  of  wavelength  for  the  same  parameters  as  in  fig.  31 
except  D  =  4000  km.   The  curves  as  in  fig.  3. 
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63 


_ 

~r  1     III 

1     1 1 

1 

1            I   .  ■ 

itit 

"T              "   n ~^ 

1            II            1       ' 

!  ■  M 

1          Mi'              '1              i    '       i       !          1 

1             '    '    1    ;            ,                  ! ,  i                 ■        / 

1  1  1 

i        1  M 

1  1 

!                    1 

,        !                                     1                1   ' 

1                                                        1 

1                            j   i 

y 

r\ 

v 

1 

1 

J' 

1 

i 

1  ,0 

! 

i          1 

^o^»- 

1 

xJt*^""^^ 

1     1 

r*4-  ''T^ "~ ^]!l" "^ ~«^«^ 

p"~-i^         "^'-.■i^lr^      " -— -a;?— ^     '    '• 

. 

''"'^*-  .^        l^^-^i-J^  1  ~      "~       _    1    _ 

!         {                                                                                                  I'M 

■^~4^      1  '^'r**^.^  i  1  1  1    ~r~"~-^;^5 

^""T-^lJ              '*""*'--»     i        '    1 

""r**^*^^  ^i_ '    I                       '      1    •  i 

1                      ""t"^    1  1  1  1  '^"*~»»-l_i 

1  1  ii~ ^^^^<ij                           1  i    '1 

0    7^       1                        i                     1                     '"l^^s^          1    '^"** 

•«L     1   i   1       i    ""i^    '                                    111:! 

1  "^^"^   !    i    ' 

1      '          '    '_.^o^                           !             II'       '   1   1 

^*"<^l    i 

1     '     i           .r'^^''^                        \                                          1     1     '           1          1 

l^' 

/i^:^S^*|          j     1                                                                         1     ' 

1                 1    ;  1  — 

r^}    j                    1    i                                    1                    1    !        '        i 

a=::S^^     i|^/ 

^sj                           I                                         1                           1         i 

iL;l;;:i=i3==^  ^^^                       L*^ 

l\      1             1                                                                ' 

f'l  f  V    '            •—  "^^^T — '                                     <^  *! 

>     i         j    1                                                                1             1 

1   lS""^^^^                                                   .^^    '    i 

:     "     \^iiX/T          it       -1-"  ^ 

Jc'>^^(tt^\  V                           /;  ZT"/  ^'J 

vXi'^  b                         1      1    1 

OS   /^<l     Mi  Or       ^^      ^r-         r\U\^}      - 

p^    ,    ^^                                                          i          1     :          i 

^■^.^^^^l        1                                  >1 

'    /            1                        ill'! 

^*^^t^^t;-A*                                              s«^ 

1  /                                1     1       i     1       1 

1        ~  ~*~                                        1^ '»»L 

\y                                                          j    j       i                  i 

1                            :    ^ — ~~^J:=;~— .— '        1     ^*k     1 

i       >^                                                              III                  1 

1                                    j                 — ~^^^^^==4=^_.^^ 

!    '~~T'V^°*' 

^^^^__'    1 

1   ix* 

''    i'"''^^^""-           1                                                                        ' 

_^,ir^ 

1       1       ^^^^T"'^                             '                                           1 ' 

n .  7  s                                                        ^>T 

_)                   I    ("N^                 1                                Mil 

•^  ,                    j    \JffiJ                 1                                ill 

.,^t--*^ '  1     '  —  •"f*'''''^ - 

^^^^jp*^                                                                      \ 

.J^ii      J^'T"^       1 

—  '■C!^*^^                                                                               i 

_,— ^"^      Jl.J'-*^         1                       1^— i^p2^ 

-  -^''^^ 

.-J^ '"ijLj^ -Til         —*— rr^i—     '        1    1 

j        1 

-,'^'''       -^- ""^  j  ^-^  =ss==^^"T'                 1 

1            1                           1 

■^T        _    ^<-«*'^_\- —                       j 

i/^.^'j^Cf-*"''''       1 

1 

yX^^ 

1    ! 

0    0^^                       1                          1                    !          1   1       1 

1    1            ,         1                               1,                             '    ■ 

15 


30 


43 


60  ir(m  sec"') 


Fig.  34.   The  profiles  of  the  zonal  \7iacl  initially  and  at  1,  2,  and 
2  1/4  days  coraputed  from  equation  (4.9)  i/it'.i  ^  =  ^^^,  q^=  0,  L  =  2000  km, 
D  =  3000  km,  B  =  30  m  sec"'  ,  A,   corresponUs  to  5  in  sec"' 
(1.59  X  10^  m^  sec"'  ),  ano  A^   =  0.   U  (0)  =  E(  l-cos(47ty/D)). 
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Fig.    35.      The   profiles   of   the   ::oaal  wind    iiiitially  and   at   2   and  2    1/4 
days   compute^!    from  equation   (4.9)  with  parameters   as    in    fig.    34   except 
A        corresponding    to      5  m  sec"'   (A^   =  A   ). 
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Fig.  36.   The  profiles  of  the  zoaal  i7ind  initially  and  at  3,  4,  and 
5  days  computed  from  equation  (4.9)  with  parameters  as  in  fig.  35 
except   D  =  6000  km  and   L  =  4000  km. 
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